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ABSTRACT 


Roll  motion  is  the  least  predictable  of  ship  motions,  largely 
because  of  viscosity  effects,  which  in  turn  causes  flow  separa¬ 
tions  to  occur.  As  a  result  of  these  phenomena,  roll  damping  is 
often  highly  nonlinear.  It  is  best  estimated  in  the  present 
state  of  the  art  by  strictly  empirical  methods. 

In  this  report,  several  relevant  flow-separation  problems 
are  studied  from  a  fundamental  point  of  view.  These  are  divided 
generally  into  two  categories:  (i)  boundary-layer  separation 
and  (ii)  sharp-edge  separation.  In  both  cases,  starting  motion 
and  oscillatory  motion  are  considered  separately.  Emphasis  is 
placed  on  rather  simple  methods  of  calculations.  Results  are 
compared  with  experiments  and/or  more  detailed  calculation 
methods  wherever  such  are  available.  Finally,  a  procedure  is 
proposed  (but  not  carried  out)  for  treating  the  formation  of 
vortex  sheets  by  a  translating  rolling  ship  hull. 
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INTRODUCTION 


The  prediction  of  hydrodynamic  forces  on  an  oscillating  ship 
hull  is  now  one  of  the  most  important  problems  of  ship  hydrodynamics. 
Among  several  factors  affecting  the  forces,  the  one  due  to  fluid 
viscosity  is  still  left  for  theoretical  treatment.  One  of  the  rea¬ 
sons  seems  to  oe  the  difficulty  of  treating  the  separating  flow 
around  a  bluff  body  like  a  ship  hull. 

There  are  several  features  of  flow  separation  around  a  ship 
hull.  In  the  absence  of  ship  speed,  the  separation  occurs  at  the 
stem  and  the  stern,  forming  vertical  vortices  during  lateral  ship 
motions.  Bilge  separation  also  occurs,  shedding  horizontal  vortices 
near  the  turn  of  the  bilge  if  the  radius  of  the  bilge  circle  is 
small  enough.  When  the  ship  moves  forward,  these  vortices  flow 
downstream,  forming  a  vortex  wake  which  is  oscillating  and  mixing 
into  the  ship  wake  caused  by  the  steady  forward  motion.  At  present, 
it  is  impossible  to  handle  all  of  these  complicated  vortex  features 
in  a  single  system.  Some  sort  of  simplification  is  necessary  to 
treat  them. 

In  this  paper,  several  attempts  are  made  to  predict  flow  sepa¬ 
ration  around  an  oscillating  ship  hull.  In  the  case  of  zero  ship 
speed,  the  problem  is  reduced  to  two-dimensional  flow  separation, 
which  can  be  treated  in  two  ways.  One  is  the  separation  from  a 
smooth  curved  surface,  which  is  discussed  in  Chapter  2.  The  other 
is  the  separation  at  a  sharp  edge  or  corner.  The  most  simplified 
case  of  the  latter,  v is.,  the  flat-plate  case,  is  treated  in  Chap¬ 
ter  3  by  a  simple  approach.  In  Chapter  4,  a  method  is  proposed  for 
obtaining  hydrodynamic  forces  on  an  oscillating  ship  hull  with 
and/or  without  forward  speed,  in  which  the  flow  is  described  in 
terms  of  vortex  sheets  shedding  from  the  separation  points  on  the 
hull.  A  numerical  calculation  scheme  for  the  vortex-sheet  forma¬ 
tion  is  stated.  However,  the  computation  is  not  completed,  partly 
because  of  the  difficulty  of  avoiding  the  self-mixing  and  the  insta¬ 
bility  of  the  vortex  sheet  in  an  oscillatory  motion,  which  is  also 
discussed  in  Chapter  4. 
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BOUNDARY-LAYER  SEPARATION 


In  order  to  consider  the  separation  problem  in  oscillatory 
motion,  it  is  easier  to  begin  the  analysis  with  a  starting  motion. 
The  motion  starts  from  rest  with  an  acceleration  in  one  direction 
and  then  ceases  with  a  deceleration  in  the  opposite  direction. 

This  might  roughly  correspond  to  an  arbitrary  single  swing  during 
an  oscillation  period. 

After  the  discussion  of  this  starting  motion  in  Section  2.1, 
the  analysis  of  oscillating  motion  according  to  Schlichting's  solu¬ 
tion  will  be  made  in  Section  2.2. 


2.1  Flow  Separation  in  Starting  Motion 


For  any  arbitrary  shape  of  body  without  sharp-edge  corners, 
the  separation  point  can  be  determined  by  boundary- layer  theory. 
Suppose  the  velocity  at  the  outer  edge  of  the  boundary  layer, 
U(x,t)  ,  is  given  in  the  form 


U  (x,  t) 


0  t  <  0  , 

u„  (x)  (t+  gt3)  ,  t  >  o  , 


) 


(2.1) 


where  x  is  measured  from  the  front  stagnation  point  along  the 
surface  of  the  body,  and  t  represents  time.  The  coefficient  8 
represents  an  additional  acceleration  (8>0)  or  deceleration  (  8 
<  0  ) .  Without  the  8  term,  the  flow  would  correspond  to  the  case 
of  a  starting  motion  with  constant  acceleration,  which  was  treated 
by  Blasius  [1] .  The  objective  of  this  section  is  to  extend  Blasius 
treatment  to  the  case  with  an  additional  deceleration  term  (  8  <  0  ) 


The  solution  can  be  obtained  by  an  expansion  in  time  series. 
Suppose  the  stream  function  ip(x,y,t)  is  expanded  in  the  form 
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=  2/vT{U50  +  +  ...  } 


Tx' 


dU, 


-  2/vt{U0t?0  +  t3  (U0BC0  +U0-^-S1  +  ...)  }  #■  (2.2) 

where  y  is  measured  normal  to  the  wall,  v  is  the  kinematic  vis¬ 
cosity,  and  the  quantities  S0'£w*  are  functions  of  a  normalized 
coordinate  n  =  y/2/vt  .  Then  the  velocity  components  u  and  v 
with  respect  to  x  and  y  are  expressed  in  the  form 


»  •  ! 

1  ■ 

V  a  - 

H  -  - 

3x‘ 


(2.3) 


The  present  governing  equation  is  the  boundary- layer  equation  for 
unsteady  motion,  which  can  be  written  in  the  form 

3  2u 


3u  .  3u  .  3u  „30  3U 
It  u3ic  +  1y  “  °7x  ~Tt 


-  0  . 

q  u  dy  - 

with  the  boundary  conditions 

u  =  v  =  0  at  n  =  0  , 
u  =  U  at  n  =  °°  . 


(2.4) 


) 


(2.5) 


Substituting  u  and  v  into  Eq.  (2.4)  and  noting  that 


= 


_ n.  ±_  f 

y=const.  n*const.  2t  3n 

we  can  obtain  the  following  time  series: 

u0<c;  +  +?»s«,(2-6) 

+  (3  (3^5  "fcS  -  j5o’-3)}  +  ...  =0  . 

The  equation  of  the  first-order  solution  S0  is  as  follows: 

5  J*  +  2  T)  5  J  ■*  4  {'g  »  —  4  #  (2*7) 

with  the  boundary  conditions 

G0(0)  -  (0)  *  0  and  ;;(»)  *  1  . 


i 


j 


-  4  - 


The  solution  was  given  by  Blasius  in  terms  of  the  velocity  profile. 


where 


1  +  ~sn  e“n  -  (l+2n2)  erfc  (n)  , 

/it 


f  n 

erfc  (n)  =  1  -  |  exp  (-n'2)  dn'. 


(2.8) 


The  second  derivative  of  ?0  at  the  wall,  which  corresponds  to  the 
skin  friction  stress,  has  the  value  ?"(0)  »  4//n  . 

The  second-order  equation  can  also  be  derived  by  rearranging 
Eq.  (2.6)  and  defining  y  =  88/uJ  : 

? +  2n  ?"  -  12  C I  -  “  4  +  4  U0'2+  S0O  +  y  (-1  +  C  o  >  /  (2.9) 

in  which  the  boundary  conditions  are  homogeneous.  Eq.  (2.9)  can  be 
solved  by  dividing  the  quantity  into  two  terms: 


?i  =  +  YSiz  * 


(2.10) 


The  solution  for  the  first  term,  ,  is  the  same  as  the  Blasius 

second-order  solution,  which  gives  the  value  ?  "  j  at  the  wall 


»  .  .  _  1  r  31  256  , 

-  7=  ll5  * 


The  equation  for  ?l2  can  be  expressed  in  the  form 


with 


?"’2  +  2  n  S"2  -  12?  J  2  =-!  +  ?; 


S12  (0)  *  ?I2  (0)  -  Cj2  (•)  -  0  . 


(2.11) 


(2.12) 


This  is  newly  obtained  here  and  its  solution  can  be  written: 


7=  e~n  (wn +r?n3 +  r^n5) 

-  erfc  (n)  [70 2  +  jn 4  +  6  » 


(2.13) 


?"2  (0)  -  3/10/F  . 


(2.14) 
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Summing  up  all  of  these  solutions,  we  obtain  the  velocity 
profile  in  the  boundary  layer: 

u  »  u0t[?j  +  +  eu;  +  8?;2) }  +  ...J  .  (2.15) 

The  skin  friction  stress,  tw  ,  is  defined  by  the  equation 

tw  =  •  (2.16) 

*  y=  0 

The  separation  point  can  be  found  by  putting  tw  »  0  in  Eq.  (2.16), 
so  that  the  condition  of  separation  is  expressed  as  follows : 

,  dU„  a 

1  +  t2  (0.427  -gJ.  +  |  3  )  =  0  .  (2.17) 

Blasius ’  original  solution  contained  only  the  first  two  terms 
in  (2.17);  the  0  term  is  new  here.  Eq.  (2.17)  can  be  interpreted 
as  determining  the  separation  time  at  the  point  where  the  velocity 
gradient  dU0/dx  and  the  deceleration  rate  0  are  prescribed. 
Another  interpretation  is  that  it  gives  the  location  of  the  separa¬ 
tion  point  where  Eq.  (2.17)  is  satisfied.  It  can  be  said  from  Eq. 
(2.17)  that  the  separation  occurs  when  the  pressure  gradient  becomes 
adverse,  and  that  the  deceleration  (  0  <  0  )  makes  the  separation 
occur  earlier  in  time  or  further  upstream  in  space,  since  a  negative 
value  of  0  corresponds  to  an  adverse  pressure  gradient  (  dU„/dx 
<  0  )  . 

We  can  apply  this  result  to  a  couple  of  cases.  First,  consider 
a  circular  cylinder  of  radius  R  oscillating  rfith  amplitude  X  and 
frequency  w  .  The  outer-edge  velocity  U(x,t)  of  the  boundary 
layer  is  given  by  potential-flow  theory: 

U  =  2Xuj  sin  ^  sin  wt  .  (2.18) 

If  we  assume  that  this  is  a  starting  motion,  then  we  obtain,  by 
equating  (2.1)  and  (2.18), 

Ufl  ■  2  Xu2  sin  ^  , 

*  -  s9u>2  . 


0 


(2.19) 
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The  deceleration  constant  6fl  can  be  obtained  by  applying  the  least- 
squares  method  over  a  swing  of  the  motion,  which  gives: 

80  =  1/8.8  .  (2.20) 

Substituting  Eq.  (2.18)  into  (2.17),  we  can  obtain  the  separation 
point  xs  at  an  arbitrary  time  instant.  For  example,  at  the  time 
t  =  7r/2u)  ,  when  the  velocity  of  the  cylinder  is  maximum,  we  obtain 
the  separation-point  angle  0S  measured  from  the  top  of  the  cylinder 

sin  0S  =  0.322/Kc  ,  (2.21) 

where  Kc  is  the  Keulegan-Carpenter  number,  7rX/R  ,  and  0S  ■ 
xs/R  -  ir/2  . 

Eq.  (2.21)  shows  that  the  separation  depends  on  the  amplitude 
of  the  motion  but  not  on  Reynolds  number  or  frequency,  according  to 
the  present  analysis.  However,  the  value  of  the  coefficient  in  Eq. 
(2.21)  is  somewhat  small  compared  with  the  experimental  value  of 
Ikeda  et  at.  [2],  as  shown  in  Fig.  2.1.  The  experimental  data  seem 
to  be  expressed  as 


sin  0 s  »  1.49/KC  .  (2.22) 

We  can  also  apply  the  present  result,  Eq.  (2.17),  to  the  case 
of  a  flat  plate  oscillating  in  its  plane,  if  we  assume  that  dUQ/dx 
=  0  everywhere  on  the  plate.  From  (2.17)  and  the  approximation 
(2.20),  we  find  that  tw  =  0  when 

wt  *  2.345  .  (2.23) 

In  this  case,  the  condition  xw  =  0  does  not  indicate  separation 
but  merely  the  turning  point  of  the  xw  oscillation.  Details  will 
be  discussed  in  the  next  section.  The  value  in  (2.23)  is  in  quite 
good  agreement  with  the  exact  solution  by  Stokes  [1] : 

3 

ojt  *  ^rr  s 


2.355  . 


(2.24) 
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This  corresponds  to  the  phase  lead,  tt/4  ,  of  the  skin  friction  with 
respect  to  the  velocity  of  the  body,  which  means  that  the  body  will 
experience  not  only  a  damping  force  but  also  an  added-mass  force  due 
to  the  fluid-viscosity  effect,  although  the  phase  is  independent  of 
viscosity. 

Through  the  present  analysis,  it  can  be  concluded  that  the  flow 
deceleration  makes  the  separation  point  move  upstream  and  it  makes 
the  separation  occur  at  earlier  time.  For  the  case  of  periodic 
oscillation,  it  corresponds  to  a  phase  lead  of  the  flow  pattern 
with  respect  to  the  velocity  of  the  body. 


2.2  Flow  Separation  in  Oscillatory  Motion 

In  the  preceding  section,  we  have  discussed  the  starting  motion 
by  expanding  the  equation  in  time  series.  Therefore  the  results 
should  be  valid  only  in  a  small  time  period,  and  its  application  to 
oscillatory  flow  is  also  restricted  within  a  small  period  of  the 
oscillation,  i.e.,  high  frequency,  or  small  displacement  of  the  body. 

It  is  interesting  to  recall  here  Schlichting ' s  approach  [1] 
to  the  periodic-oscillation  problem  of  the  boundary-layer  equa¬ 
tion.  After  following  his  procedure,  we  shall  discuss  the  separa¬ 
tion  condition. 

The  expansion  is  based  on  the  assumptions 

0  Is/lt  *  Umax/“D  “  X/D  «  1  -  (2.25) 

where  X  is  the  amplitude  of  oscillation  and  D  is  the  reference 
dimension  of  the  body.  Thus  this  is  basically  a  small-amplitude 
expansion.  Suppose  the  velocity  U(x,t)  outside  the  boundary  layer 
(the  "edge  velocity")  is  expressed  in  the  form 

U(x,t)  =  U(x)  eiwt  .  (2.26) 

Expanding  the  velocity  u  in  the  boundary  layer  in  the  series 


u  »  u0(x,y,t)  +  Uj (x,y,t)  + 


(2.27) 


-  9  - 
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putting  this  into  the  boundary-layer  equation,  (2.4),  and  applying 
the  assumption  (2.25),  we  can  obtain  the  equations  for  uq  and  ut  : 


3u0 

TT 

3ut 

TT 


-  v 


-  V 


3  zu( 

IF 

32Uj 

IF 


3U 

3t  ' 


U  ti  "  uo 


3u 
o  1y~  ' 


(2.28) 

(2.29) 


with  the  usual  boundary  conditions. 


The  solutions  of  Eqs.  (2.28)  and  (2.29)  were  given  by  Schlicht- 
ing  in  the  form 

-  (l+i)n 


c 

© 

II 

u  eiu)t{l  - 

U1  ■ 

!  Ul»  +  UM 

DU'  r  3  j 

u  l  a  = 

:  —  f'T  + 

UU'  2iu)t  r 

u  1  1  = 

- e  { 

/it  1 

1  , 


(2.30) 

(2.31) 


-2ti  __  . 

—  +  e"n  (jcos  n  +  2  sin  n) 

+  j  ®”n (sin  n  -  cos  n ) }  ,  (2.32) 

_i  e-(l.i)/In  +  +i^in)e*<l+i,n}  , 

(2.33) 


where  n  «  y/uj/2v  and  U'  =  dU/dx. 


(2.34) 


The  most  important  feature  of  these  solutions  is  that  there  appear 
(i)  a  steady  flow  in  Eq.  (2.32)  and  (ii)  a  second-harmonic  term  in 
Eq.  (2.33),  both  due  to  higher-order  effects  of  the  oscillation. 


Wow  we  have  to  pay  attention  to  the  skin  friction  variation 
on  the  wall,  which  leads  to  the  condition  of  separation.  The  skin 
friction  tw  can  be  expressed  in  the  form 


lw 


p/voj 


,i  (cot+Tr/4) 


r 

3u> 


fj+ 


(/J-lle1  (2wt-ir/4)  ^ 


(2.35) 


For  convenience  in  comparing  this  to  the  result  of  the  previous 
section,  suppose  the  edge  velocity  is  prescribed  in  the  form 


U 


U sin  wt  . 


(2.36) 


-  10 


Then  the  skin  friction  can  be  written  as  follows: 


Lw 


=  p/voj  U 


sin  (cot  +  j) 


{j-  (/2-1)  sin  (2wt+J)}  .  (2.37) 


The  first  term  coincides  with  that  of  flat-plate  oscillation,  while 
the  second  term,  which  is  always  positive,  represents  the  modifica¬ 
tion  of  the  pressure  gradient  or  of  the  amplitude  gradients  of  the 
external  flow  along  the  wall. 


Putting  tw  *  0  into  Eq.  (2.37),  we  can  derive  the  separation 
condition: 


1  dU  _  _ /?  sin  (wt  +  tt/4) _  -2  3g. 

w  3x  1/2-  (/S-l)  sin(2wt  +  ir/4) 

Eq.  (2.38)  shows  that  the  separation  occurs  when  and  where  the  value 
of  the  left-hand  side  exceeds  that  of  the  right-hand  side  in  one 
swing.  The  nature  of  the  separation  condition  is  almost  the  same 
as  in  the  previous  section  in  that  it  mainly  depends  on  the  amplitude 
of  the  motion  but  not  on  the  frequency  or  on  Reynolds  number. 

We  can  also  apply  Eq.  (2.38)  to  the  sway  motion  of  a  circular 
cylinder  of  radius  R  .  From  Eq.  (2.18)  and  Eq.  (2.36),  the  velocity 
amplitude  is  expressed  as 

U  *  2  X  id  sin  x/R  .  (2.39) 


Let  us  consider  again  the  problem  of  obtaining  the  separation  point 
xs  or  the  angle  0S  *  xg/R  -  ir/2  at  the  instant  wt  *  ir/2  .  Sub¬ 
stituting  Eq.  (2.39)  and  wt  *  tt/2  into  Eq.  (2.28),  we  obtain 

sin  9S  *  1,v98-  (  Kc  -  itX/R  )  .  (2.40) 

Eq.  (2.40)  has  the  same  form  as  Eq.  (2.21),  but  the  coefficient 
differs  much.  The  experimental  value  in  Eq.  (2.22)  is  almost  the 
average  of  Eqs.  (2.21)  and  (2.40),  as  shown  in  Fig.  2.1. 

Fig.  2-2  shows  the  comparison  of  the  separation  condition  in 
an  arbitrary  time  within  one  swing.  The  values  from  Eq.  (2.38)  and 


-  11  - 


,Eq.  (2.17) 

starting  motion  (expanded  Blasius*  method  with 
U  =  U0(t+8t3) ) 

=  U(0Jt+8rtu)3t3) 


W\ 

\\ 

\ 

\ 

\ 


Eq.  (2.38) 
oscillatory  motion 
U  =  Usinoit 


.approximation  Eq.  (2.41) 


Figure  2-2.  Condition  of  separation  in  one  swing 


and  from  an  equation  derived  from  (2.17)  have  a  similar  tendency  in 
the  middle  part  of  the  one-swing  period.  However,  at  both  ends  of 
the  period  there  is  a  difference  due  to  the  time-history  effect. 

Moreover,  an  interesting  thing  is  that  both  curves  go  to  zero 
at  wt  »  3ir/4  .  This  means  that  at  this  time  the  skin  friction 
becomes  zero  at  the  point  where  dff/dx  »  0  .  In  the  case  of  a  flat 


plate. 


3  at  the  point  where  du/dx  »  0  .  In  the  case  of  a  flat 
changes  its  sign  at  u>t  *  3rr/4  ,  while  in  the  case  of  a 


L__x 
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body  with  round  bluff  corner  the  separation  point  moves  upstream 
into  the  accelerated  flow  region  (  dU/dx  >  0  )  on  the  body  surface. 


Fig.  2-2  also  shows  an  approximation  curve  which  is  deter¬ 
mined  in  such  a  way  that  it  agrees  with  the  experiment  at  the  time 
wt  *  ir/2  ,  i.e.,  Eq.  (2.22).  The  form  of  the  curve  is  expressed 
as  follows: 


1  dU 
u  3x 


4  •  213  — j* 

loi2t2 


(2.41) 


Eq.  (2.41)  would  be  useful  to  obtain  the  separation  point  of  an 
arbitrarily  shaped  body  at  any  instant  of  one  swing  period. 


Throughout  the  preceding  analysis,  we  have  considered  that  the 
separation  is  expressed  by  the  condition  tw  =  0  .  However,  this 
condition  may  also  imply  attachment  or  a  simple  turning  of  tw 
without  separation  of  the  flow.  An  additional  condition  to  distin¬ 
guish  separation  from  the  others  can  be  found  in  the  following  way. 


Suppose  the  stream  function  i|>(x,y,t)  close  to  the  wall  is 
expanded  in  the  form 

\fj(x,y,t)  =  ay2  +  by3  +  ...  .  (2.42) 


The  coefficient  a  can  be  related  to  the  skin  friction  tw  ,  and 
the  other  coefficient,  b  ,  can  also  be  found  from  the  leading-order 
term  of  the  boundary- layer  equation,  (2.4),  expanded  near  the  wall. 
Then  Eq.  (2.42)  becomes 

'2  -  -  -3U  ,  „  3Ui 


&  -  $y  CH  ♦  o®!  • 


(2.43) 


Near  the  point  tw  *  0  ,  tw  can  be  expanded  in  the  form 


lw 


<*-xs)  +  TT  (t-ts) 


•  •  •  • 


(2.44) 


At  the  time  t  =  tg  ,  the  stream  function  ip  becomes 


(2.45) 


13 


The  quantity  in  brackets  shows  that  a  dividing  streamline  starts 
at  the  point  x  ■  xs  (  tw  »  0  ) ,  and  its  inclination,  tarn  a  ,  is 
expressed  as 


A  positive  value  of  tan  ot  means  that  the  flow  is  upwards,  i.e., 
there  is  a  separation.  A  negative  value  corresponds  to  a  reattach¬ 
ment,  and  zero  value  corresponds  to  simple  turning  of  tw  . 

For  the  separation  in  one  swing  period  (  U  >  0  )  of  oscillatory 
motion,  we  can  safely  assume  that 


0 


(2.47) 


3U 

Tt 


< 


0 


(2.48) 


Eq.  (2.48)  is  almost  satisfied  by  the  u'  value  of  Eq.  (2.38). 
Substituting  Eq.  (2.37)  into  Eq.  (2.47),  we  can  obtain  an  additional 
condition  for  separation  in  oscillatory  flow  in  the  following  simple 
form: 


0 


(2.49) 


On  the  contrary,  u"  >  0  shows  an  attachment.  Eqs.  (2.41)  and 

(2.49)  can  be  used  to  obtain  the  location  of  the  separation  point 
at  a  given  instant. 


From  the  present  analysis  of  Schlichting's  solution  for  an 
oscillatory  boundary  layer,  it  can  be  concluded  that  the  separation 
point  is  determined  by  the  derivative  of  the  edge  velocity  along  the 
surface,  i.e.,  the  pressure  gradient.  Before  the  end  of  one  swing, 
the  separation  point  moves  upstream  across  the  point  of  zero  pres¬ 
sure  gradient.  These  features  are  the  same  as  those  obtained  in 
the  preceding  section. 
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When  we  use  potential  flow  theory  to  obtain  the  edge  velocity 
distribution  along  the  surface,  we  can  roughly  predict  the  location 
of  the  separation  point  by  applying  the  criterion  formula  in  this 
section.  However,  more  precise  treatment  of  the  outer  and  inner 
flows  would  be  necessary  to  obtain  an  accurate  prediction  of  the 
separation  and  the  attachment  points.  The  present  analysis  is  a 
starting  step  for  the  problem  of  flow  separation  from  a  smooth 
surface . 


3. 


SHARP-EDGE  SEPARATION 


The  problem  of  the  separation  at  a  sharp-edge  corner  is  also 
important  for  estimating  hydrodynamic  forces  on  an  oscillating 
bluff  body  like  a  ship.  In  this  chapter#  an  attempt  is  made  to 
obtain  the  forces  in  a  quite  simple  way.  Instead  of  expressing 
the  separated  flow  in  terms  of  vortex  sheets  shed  from  the  edges, 
an  assumption  is  made  that  one  or  two  concentrated  vortices  can 
be  used.  In  the  following  two  sections,  the  motions  of  a  flat 
plate  normal  to  its  plane  are  discussed. 


3.1  Starting  Motion  of  a  Flat  Plate 

To  obtain  the  normal  force  on  a  flat  plate  starting  from  rest 
with  a  constant  velocity,  the  following  simple  procedures  are 
adopted : 

i)  Put  a  pair  of  vortices  at  the  midpoints  of  the  path  lines 
of  the  separating  edges  (see  Fig.  3.1). 

ii)  Determine  the  vortex  strength  by  Kutta  condition. 


To  begin,  suppose  that  the  flat  plate  of  breadth  2  in  the 
z  plane  is  mapped  into  a  segment  of  length  2  in  the  5  plane, 
as  shown  in  Fig.  3.1.  The  edge  z  =  i  corresponds  to  5  =  0. 

The  mapping  function  is 

z  *  /jf2  - T  .  (3.1) 

Then  assume  that  a  vortex  of  strength  r  lies  on  each  of  the  path 
lines  of  the  separating  edges  at  an  arbitrary  distance  x0  from 
the  edge: 


or 


+  i 


;  -  ;0  -  50  +  in0 


(3.2) 
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Figure  3.1.  Assumed  pair  of  vorticies  behind  a  flat 
plate  starting  from  rest 


in  the  £  plane.  The  two  vortices  have  opposite  directions  of 
rotation. 


jugate  velocity,  that  is, 

q(z)  *  u (x,y)  -  iv(x,y)  *  QU)/j|  , 

where 

Q<c>  -  v  +  ~  • 


i  to  the  plate. 

there 

i  shown  in  Fig. 

3.1. 

of  the  [complex]  con- 

QU)/§§  , 

(3.3) 

(3.4) 

At  the  edge  point  z  =  i  or  c  =  0  ,  the  Kutta  condition  is 
applied  to  determine  the  vortex  strength  r  : 


r  =  -  2ttV  1  ioi° 
'•So  " 

✓5T 

*  -  2ttV — /x?  • 


-  2irV272  /Xo  +  4  ^X°  +  /xo  +  I//2 


For  small  values  of  x0  ,  this  becomes  approximately 

r  a  -  2  TlV  /x7  . 


(3.5) 

(3.6) 

(3.7) 


The  vortex  is  assumed  to  lie  on  the  midpoint  of  the  path  length  s 
of  the  edge: 


s  1 
7  7 


v  t  . 


(3.8) 


The  force  on  the  plate  is  easily  derived  from  the  impulsive 
momentum  I  of  the  flow  field, 

F  -  a|  =  ^{p*b2V-2py0r}  ,  (3.9) 

where  b  =  half-breadth  =  1  ,  y0  =1  .  Then 

F  -  -  2p|I  .  (3. 10) 

From  Eqs.  (3.7)  and  (3.8),  the  first  approximation  of  r  for  small 
x0  is 

r  *  -  2  TrV/\7t7T  .  (3.11) 


The  force  can  be  calculated  in  the  form 


C 


N 


F 

(1/2) p • 2bV* 


/Itt/ZF  , 


(3.12) 


with  t'  =  Vt  .  Figs.  3.2  and  3.3  show  results  calculated  from 
Eqs.  (3.11)  and  (3.12).  Eq.  (3.11)  becomes  quite  close  to  Wede- 
meyer's  exact  solution  [3].  Eq.  (3.12)  is  also  acceptable  compared 
to  Kudo's  calculation  [4], 


-r 

2itVb 


O:  /t"'/2  ,  Eq.  (3.11) 
A:  Eq.  (3.6) 


Figure  3-2.  Vortex  generation  in  starting  motion 
of  flat  plate  of  half -breadth  b*l. 


o 


Figure  3-3.  Normal  force  on  flat  plate  m 
starting  motion 

The  success  of  this  simple  approach  can  be  understood  in  the 
following  way:  The  velocity  qv  at  a  vortex  position  can  be 
obtained  by  dropping  the  term  corresponding  to  that  vortex.  Thus, 


qv  =  Qv/|f 


(3.13) 


vfl  +  ^0  = 

(Cn-C‘n)2J 


x0/x^  T 4  j 

2  (x0  /xTTT  -  x? )  J 


0  0 


o  r  o 

■oKXo 


So  we  find  that 


(3.14) 


(3.15) 


Eq.  (3.15)  means  that  the  vortex  velocity  is  half  of  the  velocity 
in  the  absence  of  that  vortex  pair.  Although  this  does  not  mean 
that  the  location  of  the  vortex  is  on  the  half-way  point  of  the  path 
of  the  edge,  this  fact  is  quite  interesting. 
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If  we  again  assume  that  the  vortex  lies  on  the  line  y0  =  1 
and  obtain  the  distance  x0  from  the  edge  by  solving  the  equation 

dx.  _ 

=  uv(x0)  =  Re{  qv} 


=  V/4/x^  for  small  x0  , 

the  result  is 


Applying  this  to  the  case  of  starting  motion. 


(3.16) 


(3.17) 


*o  = 


(3.18) 


and  substituting  this  into  Eq.  (3.7),  we  obtain 

)  1/3 

Vt  .  (3.19) 

Eq.  (3.19)  is  also  shown  in  Fig.  3.2  and  gives  slightly  lower  values 
than  Eq.  (3.11)  in  the  range  Vt  >  1  . 


r  =  -  2  -a  v 
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On  the  other  hand,  Eq.  (3.6),  the  original  expression,  gives 
higher  values  when  x0  is  taken  as  s/2  ,  the  midpoint.  Therefore 
we  can  realize  that  Eq.  (3.11)  compensates  two  errors,  one  arising 
from  the  assumption  of  vortex  location  (  x0  =  s/2  ) ,  and  the  other 
coming  from  the  simplification  of  the  r  expression. 

There  seems  to  be  no  evidence  for  the  present  assumption,  Eq. 
(3.8),  on  the  vortex  location.  However,  we  can  cite  here  some  sup¬ 
port,  for  instance,  Pierce's  photograph  [5],  Bollay's  assumption 
that  the  trailing  vortex  of  a  low-aspect-ratio  wing  sheds  at  half 
the  angle  of  incidence,  and  the  numerical  calculations  made  by  Fink 
and  Soh  [ 3 ]  and  Kudo  [ 4 ] . 

In  conclusion,  it  can  be  said  that  the  present  simple  approach 
predicts  the  force  on  a  flat  plate  with  good  accuracy.  However,  it 
can  not  be  expected  that  this  method  will  also  give  a  valid  detailed 
description  of  the  flow  field. 
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3.2  Oscillatory  Motion  of  a  Flat  Plate 

According  to  the  observations  by  Ikeda  et  al  [2]  and  the 
numerical  calculations  by  Kudo  (4]  and  by  Fink  and  Soh  [5],  the 
flow  around  a  flat  plate  oscillating  normal  to  its  plane  is  quite 
interesting. 

Roughly  speaking,  a  pair  of  vortex-sheet  cores  is  created 
during  each  swing  of  the  body.  However,  toward  the  end  of  a  swing, 
before  the  plate  motion  ceases ,  the  sign  of  the  vorticity  density 
near  the  edge  becomes  opposite  to  that  of  the  vortex  core  created 
during  the  swing.  This  means  that  the  strength  of  the  shed  vor¬ 
ticity  has  a  phase  lead  with  respect  to  the  velocity  of  the  plate. 

The  physical  interpretation  of  this  fact  seems  to  be  that  the 
velocity  at  the  rear  side  of  the  edge  induced  by  the  vortex  core 
just  created  becomes  greater  than  the  velocity  on  the  front  side 
near  the  edge. 

During  the  next  swing,  the  previous  vortex  moves  downstream 
past  the  plate  and  a  new  vortex  core  is  generated  behind  the  plate. 
According  to  the  calculations  by  Kudo  or  by  Fink  and  Soh,  the  loca¬ 
tion  of  the  previous  vortex  is  nearly  above  the  newly-created  vortex. 
During  continuous  oscillatory  sway  motion,  the  vortices  thus  created 
seem  to  flow  away  from  the  plate. 

In  order  to  formulate  this  kind  of  flow  mathematically,  we 
can  again  apply  the  idea  in  the  preceding  section.  With  the  coor¬ 
dinates  fixed  on  the  body,  the  upstream  velocity  can  be  expressed 

V  *  V  sin  wt  ,  V  *  u)X  .  (3.20) 

Assume  that  the  flow  field  is  represented  by  two  pairs  of  concen¬ 
trated  vortices  which  lie  somewhere  on  the  horizontal  lines  through 
the  edges.  Suppose  their  distances  from  the  edge  are  x0  and  x, 
and  their  strengths  r0  and  Tj  ,  as  shown  in  Fig.  3.4.  When  the 
flow  goes  right  (  V  >  0  ) ,  we  assume  that  r„  represents  the  pre- 
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vious  vortex,  which  is  constant, 
and  Tj  the  new  one,  which  is 
growing  with  time. 

Assume  that  the  total  circu¬ 
lation  in  the  upper  half  plane  is 
of  the  form 

T  =  T  cos  (tot  +  £ )  ,  (3.21) 

with  the  amplitude  r  and  the 
phase  lead  e  .  At  the  time  tot 
=  -  e  ,  there  is  only  a  pair  of  vo: 
tices,  the  strength  of  which  is 

r0  =  r  •  (3.22)  Figure  3.4.  Two-vortex-pair 

model  for  oscillating 

Thereafter  F0  is  constant.  How-  flat  plate 

ever,  a  new  pair,  of  strength  , 

grows  with  time  until  it  reaches  its  maximum  value, 

Tj  =  -  2  T  at  tot  =  it  -  e  .  (3.23) 

At  cot  =  tt  -  e  ,  the  total  circulation  T  equals  -T  ,  that  is, 

T  =  Tg+Tj  =  -  T  at  tot  =  tt  -  £  .  (3.24) 

However,  the  flow  field  should  again  be  expressed  by  one  pair  of 
vortices  of  strength  Tj  at  cot  =  tt-e,  so  that  the  previous  vor¬ 
tex  should  be  mixed  into  the  new  vortex  r  at  tot  =  it  -  e  .  In 
other  words,  T0  has  traveled  along  an  unknown  route  and  joined 
with  the  new  vortex  r  l  at  <ot  =  it  -  e  .  In  this  model ,  we  do  not 
concern  ourselves  with  the  path  of  the  vortex  T0  .  Instead,  we 
have  assumed  the  variation  of  F  as  in  Eq.  (3.21). 

To  determine  the  vortex  strengths,  we  apply  the  Kutta  condition 
at  the  edges  of  the  plate.  At  cot  *  -  e  ,  there  is  a  vortex  pair 
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r  0  .  The  velocity  at  the  edge  z  *  i  or  5=0  can  be  expressed 
in  the  5  plane  by  using  the  same  mapping  function  as  in  Eq.  (3.1): 


Q  = 


v  + 


(3.25) 


5  =  0 

The  Kutta  condition (  Q  =  0  )  requires  that 


r  =  -  2iTiv  . 

5o  -  50 

For  small  values  of  x0  ,  this  becomes 

r0  =  27T/|x0T  V  sin  e  =  T  . 


(3.26) 


(3.27) 


At  times  such  that  wt  >  - e  ,  there  are  two  vortex  pairs ,  F 
and  T  x  .  The  Kutta  condition  for  the  growing  vortex  F l  can  be 
written  in  the  form 


or 


Q 


v (t)  +  v 0  (t)  + 


-  0  , 


rx  =  -  2,rrv/x^"[V+ V0  ]  for  small  xx  , 


(3.28) 

(3.29) 


where  VQ  represents  the  induced  velocity  of  the  previous  vortex, 
the  strength  of  which  (  ro  )  is  now  constant.  But  the  location  of 
r0  is  now  unknown,  so  that  VQ  is  unknown. 

For  the  growing  vortex  r  ,  we  assume  the  velocity  xx  as 
follows : 

Xj  =  y[V+V„]  for  -e  <  wt  <  tt  -  e  .  (3.30) 

The  quantity  V +  V„  is  the  velocity  at  5  =0  in  the  5  plane, 
which  means  that  the  incoming  velocity  at  the  edge  has  increased  by 
the  amount  VQ  due  to  the  effect  of  the  previous  vortex.  This 
assumption  is  similar  to  the  one  that  the  vortex  lies  on  the  mid¬ 
point  of  the  travel  path  of  tne  edge. 


From  Eqs.  (3.28),  (3.29),  and 

r,  =  r-ro  =  7{cos  ut+  o  -  1}  ,  (3.31) 
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we  can  obtain  xL  : 

yirx|/2  =  £{u>t  +  e  -  sin  (uit  +  £)} 


(3.32) 


It  should  be  noted  that  we  have  not  assumed  the  form  of  VQ 
the  instant  wt  =  ir  -  e  ,  Eq.  (3.32)  becomes 


At 

(3.33) 


and,  from  symmetry  of  the  vortex  locations  between  the  times  tut  = 
-e  and  ir  -  e  ,  we  obtain  another  relation  from  Eq.  (3.27): 


x. 


l2irV  sin  eJ 

Eqs.  (3.33)  and  (3.34)  can  be  arranged  into  the  forms: 

1/2 


r  = 


2n  sin  e 
3ir 


'  3/2  f3] 

18J 


; 


•  (sin  e)  X 


(3.34) 

(3.35) 

(3.36) 


Note  that  Xj  or  x0  is  now  the  vortex  location  at  the  time  when 
its  strength  becomes  maximum,  that  is,  at  ait  =  e  or  ir  -  e  . 


Assumption  on  the  phase.  The  phase  lead  e  of  the  circulation 
r  with  respect  to  the  velocity  V  plays  an  important  role  in  this 
analysis.  It  is  introduced  to  express  an  effect  of  the  previous 
swing.  However,  it  is  still  not  clear  what  is  a  reasonable  condi¬ 
tion  to  determine  the  value  of  £  .  Here  we  assume  the  value  to  be 


£  *  J  =  45°  .  (3.37) 

The  support  of  this  assumption  can  be  found  in  several  examples  of 
viscous  flow  theory.  For  instance,  the  shear  stress  in  the  oscilla¬ 
tory  boundary  layer  has  the  same  phase.  Kudo's  numerical  calculation 
also  showed  a  phase  lead  of  about  it/4  in  the  case  of  flat-plate 
sway.  We  can  expect  that  the  phase  lead  ir/4  might  be  derived  from 
a  high-frequency  expansion  of  the  equation  system  for  this  kind  of 
flow.  The  problem  is  left  for  future  study. 
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The  force  can  be  obtained  by  taking  the  time  derivative  of 
the  impulsive  momentum  of  the  flow  field: 


F 


dV 

p1T3t 


-  2 


n  dF 

p3t 


(3.38) 


Thus, 


F  =  p{irXu>2  +  2BXw2*/x  sin  e}cos  out 


+  2pf3Xu)2  /£  cos  e  sin  ut  , 


(3.39) 


where 


) 3/2  f  3) 

3  =  2  it  sin  e 


3)1/2 


=  5.735  for 


e  =  tt/4  . 

(3.40) 


In  this  analysis,  the  force  has  only  simple-harmonic  components 
at  the  fundamental  frequency  of  oscillation.  In  order  to  compare 
tiie  result  with  other  results,  we  have  to  transform  the  damping 
term  into  the  usual  nonlinear  form: 


F  =  Cmp7rXa)2  cos  ut  +  CD*^-*2  •X2u)2sin2a)t  .  (3.41) 

Equating  the  work  done  in  one  swing  period,  as  derived  from  Eqs. 
(3.39)  and  (3.41),  we  obtain 

CD  =  jTrecos  £  =  9v|55-  for  e  =  tt/4  ,  (3.42) 

Cm  *  1  +  sin  e  /X  =  1  +  2.582/x  for  e  *  tt/4  .  (3.43) 

Figures  3.5  and  3.6  show  the  comparison  of  Eqs.  (3.42)  and  (3.43) 
with  the  numerical  calculations  by  Kudo,  and  the  Keulegan  and 
Carpenter  experiments.  The  agreement  of  the  CD  values  is  quite 
good,  while  Cm  agrees  with  Kudo's  calculation,  which  is  about 
twice  as  large  as  the  value  from  the  Keulegan  and  Carpenter  experi¬ 
ments.  The  reason  for  this  discrepancy  is  still  unknown. 


The  conclusion  of  this  section  can  be  summarized  as  follows: 

A  simple  approach  to  the  problem  of  a  flat  plate  oscillating  normal 
to  its  plane  is  attempted.  The  flow  field  is  represented  by  two 
pairs  of  vortices.  The  variation  of  the  vortex  strength  is  assumed 
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to  be  sinusoidal  and  the  vortex  velocity  to  be  half  the  flow 
velocity  at  the  separating  edge.  An  additional  assumption  on  the 
phase  lead  of  the  vortex  strength  is  also  made.  With  the  Kutta 
condition  satisfied  at  the  edge,  the  results  are  expressed  in 
simple  forms  and  agree  with  other  works,  especially  for  Cq  . 


Figure  3.5.  Added-mass  coefficient  of  oscillating  flat  plate 


1  2  3X4  5 


Figure  3.6.  Drag  coefficient  of  oscillating  flat  plate 


4.  VORTEX-SHEET  FORMATION  AROUND  SHIP  HULL 

This  chapter  describes  an  outline  of  an  idea  and  a  procedure 
for  the  numerical  computation  of  vortex  formation  around  an  oscil¬ 
lating  ship  hull  with  and  without  forward  speed.  The  method  is 
based  on  the  slender-body  assumption  which  enables  us  to  treat  the 
problem  in  the  two-dimensional  plane  of  the  cross-sections. 

However,  a  difficulty  has  arisen  in  the  numerical  computation, 
not  because  of  a  defect  of  the  basic  concept  but  mainly  because  of 
the  difficulty  of  treatment  of  the  self-crossing  of  the  vortex 
sheet  in  an  oscillatory  motion.  This  situation  is  also  described. 

Although  the  numerical  results  have  not  been  obtained,  it  is 
worthwhile  to  state  the  procedure  briefly. 


4.1  Basic  Assumption  and  Procedure 

The  objective  is  to  obtain  the  flow  field  around  a  ship,  as 
well  as  the  force  and  moment  when  the  ship  is  moving  forward  and 
oscillating  in,  say,  roll.  Neglecting  wave  effects  and  skin  fric¬ 
tion,  we  reduce  the  problem  to  that  of  predicting  the  separating 
flow  around  the  hull.  This  can  be  represented  by  a  vortex-sheet 
system  shedding  from  separation  points  on  the  hull  if  we  neglect 
the  viscous  diffusion  of  vorticity.  Recent  calculations  on  the 
vortex-sheet  configuration  in  unsteady  two-dimensional  flow,  carried 
out  by  Fink  and  Soh  [3]  and  Kudo  [4],  showed  us  that  this  kind  of 
treatment  can  predict  the  forces  on  the  body  fairly  well.  Thus  we 
can  expect  that  this  approach  could  be  applied  to  a  ship  hull  form 
providing  that  the  separation-point  condition  is  treated  in  an 
appropriate  way. 

The  following  three  basic  assumptions  are  made  in  the  present 
method: 


26  - 


i)  Only  longitudinal  vortices  are  shed  at  the  shedding  points 
in  the  cross-sectional  plane  (see  the  coordinate  system  in  Fig. 
4.1).  According  to  Fuwa's  [6]  analysis,  which  is  based  on  the 
slender-body  assumption,  the  problem  can  be  treated  in  the  two- 
dimensional  plane  of  each  section.  In  this  case,  the  longitudinal 
component  wx  of  the  vorticity  is  dominant,  so  that  <jy  ,  co2  « 
o(wx)  and  y  ,  2  =  o(x)  . 

ii)  The  shed  vortex  is  convected  downstream  with  the  forward 
velocity  U  without  change  of  strength.  This  assumption  might 
seem  to  be  rough  compared  to  that  of  the  slender-body  theory.  It 
should  be  considered  as  a  first  step  for  the  vortex  convection  in 
the  longitudinal  direction. 


iii)  The  separation  point  is  treated  in  two  ways,  as  a  sharp- 
edge  separation  and  as  a  smooth-surface  separation.  The  former 
can  be  treated  in  the  same  way  as  in  the  recent  works  [3]  [4].  For 
the  latter,  the  position  of  the  separation  point  can  be  predicted 
by  boundary- layer  theory,  with  the  potential-flow  velocity  field 
as  described  in  a  previous  chapter  of  this  paper. 


x 
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The  general  procedure  is  to  carry  out  a  time  integration  of 
the  vortex-path  equation  by  using  the  velocity  field  at  the  instant, 
and  to  find  the  location  and  the  strength  of  a  newly  shed  vortex  at 
each  time  step.  To  begin  with,  the  following  three  pre-calculations 
are  needed: 

i)  Representation  of  the  hull  form  by  a  mapping  function 
transforming  the  section  shape  into  a  circle. 

ii)  Calculation  of  the  potential-flow  field  around  the  hull 
in  the  presence  of  ship  forward  speed  and  roll  velocity,  using 
slender-body  theory. 

iii)  Calculation  of  the  separation  point  in  each  cross-section 
as  a  function  of  time  t  in  one  roll  swing,  using  the  girthwise 
distribution  of  crossflow  velocity  amplitude  and  its  derivatives. 

When  the  value  of  the  first  derivative  is  large  enough,  the  point 
is  taken  as  a  sharp  edge. 

The  time-step  integration  is  carried  out  in  the  following  way: 

iv)  Calculation  at  the  shedding  point  to  determine  the  velocity, 
the  direction,  and  the  strength  of  the  newly-shed  vortex  sheet. 

v)  Calculation  of  the  vortex-sheet  velocity. 

vi)  Calculation  of  the  moment  acting  on  the  ship. 

vii)  Time  integration  to  obtain  the  new  position  and  strength 
of  the  vortex  sheet  at  the  time  after  one  time  step.  The  procedure 
from  iv)  to  vii)  is  repeated  for  several  swing  periods  of  the 
motion. 


4 . 2  Formulation 

With  the  coordinates  as  in  Fig.  4.1,  the  ship  section  shape 
can  be  represented  by  a  Lewis  form: 
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x  *  y  +  iz 

-  a,;  +  a2;"1  +  a 3;~3  .  (4.1) 

The  velocity  potential  $(x,y,z,t)  is  expressed  in  the  form 

$  *  Ux  +  i?e  {UF ,  +  ftF2  +  Fj]  ,  (4.2) 

where  U  and  ft  are  the  ship  forward  velocity  and  the  roll 
angular  velocity,  respectively.  The  complex  potential  F,  at 
each  section  corresponds  to  the  forward  motion  and  is  obtained  from 
Tuck  and  von  Kerc2ek's  expression  [7].  The  potential  F2  due  to 
the  ship  roll  motion  is  also  obtained  by  assuming  that  the  double 
body  of  the  section  rotates  around  the  point  0  .  This  can  be 
written 


F2  =  -i{a2 (a,  +  a3) ;”2  +  a1a35”1*}  . 


(4.3) 


The  potential  F3  represents  the  vortex-sheet  potential  which  is 
expressed  in  the  form 


i  fs 

5?r  j 


2 

log 


Ci  -  C  +  Ci 


!i 


5-1/;,  ;  +  !/;, 


y(Cx)  ds  , 


(4.4) 


where  lies  on  the  vortex  sheet  where  the  vorticity  density  is 

y  ,  and  ds  is  the  length  element  of  the  sheet.  The  integration 
should  cover  all  sets  of  vortex  sheets  in  the  double-body  section. 


The  shedding  condition  can  be  given  in  the  following  ways: 
When  separation  occurs  on  the  smooth  surface,  the  location  of  the 
separation  point  can  be  derived  from  the  boundary-layer  theory, 
i.e. ,  from  Eq.  (2.41),  at  a  prescribed  time  t  during  the  swing. 
The  vortex  sheet  is  assumed  to  shed  tangentially  to  the  surface, 
so  that  the  density  y  and  the  velocity  qvs  of  the  vortex  sheet 
can  be  given  as 


Y 


ISrtl 


<3vs 


1  _ 
qrt  ' 


(4.5) 


where  qrt  represents  the  relative  tangential  velocity  at  the 
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separation  point,  which  can  be  obtained  approximately  from  the 
potential  F2  .  When  the  separating  edge  is  sharp,  the  shedding 
condition  is  given  in  a  manner  similar  to  that  of  Fink  and  Soh's 
treatment,  in  which  the  density  is  determined  by  a  Kutta  condition 
at  the  edge  and  the  velocity  is  calculated  at  a  point  a  small 
distance  away  from  the  edge. 

Since  it  is  not  easy  to  transform  an  arbitrary  point  in  the 
real  plane  into  the  mapped  5  plane,  we  carry  out  all  of  these 
calculations  in  the  5  plane,  noting  that  in  the  transformation 
the  strength  of  the  circulation,  yds  ,  remains  constant,  so  that 
the  density  changes  by  the  rate  of  |dx/dc|  . 

In  the  numerical  calculation,  the  vortex  sheet  is  discretized 
by  a  finite  number  of  vortex-aheet  segments  of  small  length,  with 
linear  variation  of  vorticity  density  on  each  segment.  The  time 
integration  is  carried  out  by  Euler's  scheme. 


4.3  Discussion 

According  to  the  procedure  briefly  mentioned  above,  the 
numerical  calculation  was  carried  out  for  the  case  of  a  flat  plate 
(  LxB*2d  “  2x0x0. 4  )  starting  to  roll  around  its  central  axis  with¬ 
out  forward  velocity.  Fig.  4.2  shows  an  example  of  the  result  of 
the  vortex-sheet  generation  relative  to  the  body  axis.  As  time 
proceeds,  a  vortex  sheet  is  created  on  the  right-hand  side  of  the 
plate  during  the  first  swing  and  then  it  flows  to  the  other  side  in 
the  second  swing.  At  a  particular  moment,  the  vortex  sheet  crosses 
itself,  and  the  calculated  moment  shows  scattering,  as  in  Fig.  4.2. 
Since  this  kind  of  vortex-sheet  crossing  is  not  expected  in  the 
real  flow,  there  must  be  a  defect  in  the  numerical  procedure. 
Several  attempts  have  been  made  concerning  the  time-step  interval, 
the  shedding-point  condition,  etc.  However,  the  difficulty  could 
not  be  removed. 


Figure  4.2.  Vortex- sheet  formation  and 
moment  variation  of  rolling  flat  plate  with  amplitude  10° 
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The  same  phenomenon  appeared  in  the  calculations  by  Kudo  [4] 
and  Fink  and  Soh  [3],  who  discretized  the  vortex  sheet  into  a 
number  of  concentrated  point  vortices.  In  those  treatments,  the 
vortex-sheet  crossing  might  not  greatly  affect  the  calculation  of 
the  force  or  the  time  integration  of  the  sheet  configuration.  In 
the  present  case,  however,  the  crossing  becomes  a  definite  diffi¬ 
culty  so  that  further  time  steps  cannot  be  achieved. 

Although  the  present  trial  calculation  has  not  been  successful, 
the  basic  idea  of  the  present  method  still  seems  to  be  worth  pur¬ 
suing.  The  performance  of  this  calculation  must  be  left  for  a 
future  study. 


! 


5. 


CONCLUSION 


Three  different  approaches  have  been  made  to  the  problem  of 
flow  separation  in  unsteady  motion  of  a  body.  In  Chapter  2,  sepa¬ 
ration  from  a  smooth  curved  surface  is  discussed  from  the  view  of 
boundary- layer  theory.  It  is  found  through  the  analyses  of  a 
starting  motion  and  of  an  oscillatory  motion  that  the  separation 
point  depends  mainly  on  the  amplitude  ratio,  providing  that  the 
ratio  is  small. 

The  problem  of  the  sharp-edge  separation  is  discussed  in  a 
quite  simple  way  through  the  use  of  a  small  number  of  vortex 
pairs.  This  method  is  applied  to  the  starting  and  oscillatory 
motions  of  a  flat  plate ,  the  motion  being  normal  to  the  plane  of 
the  plate.  The  force  is  expressed  in  a  simple  formula,  which 
agrees  well  with  experiments  and  with  others'  numerical  computa¬ 
tions  . 


In  the  last  chapter,  a  procedure  to  obtain  the  vortex-sheet 
formation  around  a  ship  hull  oscillating  and  moving  forward  is 
described.  Although  the  numerical  calculation  has  not  been  com¬ 
pleted,  because  of  the  difficulty  of  the  numerical  treatment  of  a 
two-dimensional  vortex-sheet  configuration,  the  procedure  itself 
seems  to  be  adequate.  A  slight  improvement  of  the  vortex  calcula¬ 
tion  might  make  the  method  work  well. 

The  author  believes  that  these  discussions  and  analyses 
constitute  a  useful  step  toward  the  final  goal  of  the  prediction 
of  ship  viscous  roll  damping. 
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1  Mare  Island  Naval  Shipyard 

Shipyard  Technical  Library  (202.3) 
Vallejo,  CA  94592 

1  Assistant  Chief  Design  Engineer 

for  Naval  Architecture  (Code  250) 
Mare  Island  Naval  Shipyard 
Vallejo,  CA  94592 

2  U.S.  Naval  Academy 
Annapolis,  MD  21402 
Attn :  Technical  Library 

Dr.  S.A.  Elder 

1  Naval  Postgraduate  School 
Monterey,  CA  93940 
Attn:  Library  (2124) 

1  Study  Center 

National  Maritime  Research  Center 
U.S.  Merchant  Marine  Academy 
Kings  Point,  LI,  NY  11024 

1  The  Pennsylvania  State  University 
Applied  Research  Laboratory  (Library) 
P.0.  Box  30 

State  College,  PA  16801 

1  Dr.  B.  Parkin,  Director 

Garfield  Thomas  Water  Tunnel 
Applied  Research  Laboratory 
P.O.  Box  30 
State  College,  PA 

1  Bolt,  Beranek  &  Newman  (Library) 

50  Moulton  Street 
Cambridge,  MA  02138 

1  Bethlehem  Steel  Corporation 
25  Broadway 
New  York,  NY  10004 
Attn:  Library  -  Shipbuilding 


1  Cambridge  Acoustical  Associates,  Inc. 
54  Rindge  Ave  Extension 
Cambridge,  MA  02140 

1  R  &  D  Manager 

Electric  Boat  Division 
General  Dynamics  Corporation 
Groton,  CT  06340 

1  Gibbs  &  Cox,  Inc.  (Tech  Info  Control) 
21  West  Street 
New  York,  NY  10006 

1  Department  of  Ocean  Engineering 
University  of  Hawaii  (Library) 

2565  The  Mall 

Honolulu,  HI  96822 

2  Institute  of  Hydraulic  Research 
The  University  of  Iowa 

Iowa  City,  IA  52240 
Attn :  Library 

Dr.  L.  Landweber 

3  Department  of  Ocean  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge ,  MA  02139 

Attn:  Prof.  P.  Leehey 
Prof.  J.  Newman 
Prof.  J.  Kerwin 

1  Engineering  Technical  Reports 
Room  10-500 

Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 

2  St.  Anthony  Falls  Hydraulic  Laboratory 
University  of  Minnesota 
Mississippi  River  at  3rd  Avenue  S.E. 
Minneapolis,  MN  55414 

Attn:  Dr.  Roger  Arndt 
Library 

2  Department  of  Naval  Architecture 

and  Marine  Engineering  -  North  Campus 
University  of  Michigan 
Ann  Arbor,  MI  48109 
Attn:  Library 

Dr.  T.  Francis  Ogilvie 
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2  Davidson  Laboratory 

Stevens  Institute  of  Technology 
711  Hudson  Street 
Hoboken,  NJ  07030 
Attn :  Library 

Dr.  J.  Breslin 

2  Stanford  University 
Stanford,  CA  94305 
Attn :  Engineering  Library 
Dr.  R.  Street 

1  Webb  Institute  of  Naval  Architecture 
Attn :  Library 
Crescent  Beach  Road 
Glen  Cove,  LI,  NY  11542 

1  National  Science  Foundation 
Engineering  Division  Library 
1800  G.  Street  NW 
Washington,  DC  20550 

1  Dr.  Douglas  E.  Humphreys,  Code  794 
Naval  Coastal  Systems  Laboratory 
Panama  City,  FL  32401 
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